EPJ manuscript No. 

(will be inserted by the editor) 



Isovector deformation and its link to the neutron shell closure 
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Abstract. DWBA analysis of the inelastic 30 ~ 40 S(p,p') and 18 ~ 22 0(p,p') scattering data measured in the 
inverse kinematics has been performed to determine the isoscalar (So) and isovector (Si) deformation lengths 
of the 2J" excitations in the Sulfur and Oxygen isotopes using a compact folding approach. A systematic 
iV-dependence of So and Si has been established which shows a link between Si and the neutron-shell 
closure. Strong isovector deformations were found in several cases, e.g., the 2+ state in 20 O where <5i is 
nearly three times larger than So- These results confirm the relation Si > So anticipated from the core 
polarization by the valence neutrons in the open-shell (neutron rich) nuclei. The effect of neutron shell 
closure at N = 14 or 16 has been discussed based on the folding model analysis of the inelastic 22 0+p 
scattering data at 46.6 MeV/u measured recently at GANIL. 
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The neutron and proton contributions to the structure 
of the lowest 2 + excited states are known to be quite dif- 
ferent in the neutron-rich nuclei due, in particular, to a 
strong polarization of the core by valence neutrons pQ . In 
the distorted-wave Born approximation (DWBA) for the 
inelastic hadron scattering, the different neutron and pro- 
ton contributions to the nuclear excitation are explicitly 
determined by the isospin dependence of the inelastic form 
factor (FF). 

In general, the isospin-dependent nucleon optical po- 
tential (OP) can be written in terms of the isoscalar (IS) 
and isovector (IV) components 2. as 

U(R) = U {R) ± eUx{R), e = (N - Z)/A, (1) 

where the + sign pertains to incident neutron and - sign to 
incident proton. While the strength of the Lane potential 
U\ has been studied since a long time in the (jp,p) and 
(n,n) elastic scattering and (p,n) reaction studies, very 
few attempts were made to study the isospin dependence 
of the inelastic FF. Within a collective-model prescription, 
the inelastic FF for the nucleon-nucleus scattering is ob- 
tained by "deforming" the OP ([T]) with scaling factors 8 
known as the nuclear deformation lengths 



i) (p, n) reaction leading to the excited isobar analog 
state. It was shown [3], however, that the two-step mech- 
anism usually dominates this process and the calculated 
DWBA cross sections are not sensitive to <5i . 

ii) Another way is to extract 5q and 5i from the (p,p') 
and (n, n') inelastic scattering measured at the same inci- 
dent energy and exciting the same target state [3l|4] . Such 
double measurements are presently not feasible with the 
beams of unstable nuclei. 

We have recently suggested a compact folding method 
[5] to determine <Wi) based on the DWBA analysis of the 
ip,p') data only. In this approach, instead of deforming 
the OP, we build up the proton and neutron transition 
densities of a 2 A -pole excitation (A > 2) by using the 
Bohr-Mottelson prescription [5] separately for protons and 
neutrons 
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Here Pg B (r) are the proton and neutron ground-state (g.s.) 
densities and S T are the corresponding deformation lengths. 
Given the explicit proton and neutron transition densities, 
one can obtain from the folding model [7] the inelastic 
proton-nucleus FF in terms of the IS and IV parts as 



The explicit knowledge of the isoscalar (<5 ) and, especially, 
isovector (Si) deformation lengths would give us vital in- 
formation on the structure of the nuclear excitation under 
study. There are only two types of experiment that might 
allow one to determine the IV deformation length 5\ based 
on the prescription ((5J: 



F(R) = F {R)-sF 1 (R), 



(4) 



where Fo(R) and Fi(R) are determined 'T\ from the sum 
(px + P\) an( i difference (p 1 ^ — p 1 ^) of the neutron and 
proton transition densities ([3]), respectively. One can see 
that Fi(R) is just the prototype of the Lane potential in 
the inelastic nucleon scattering. It is natural to represent 
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the IS and IV parts of the nuclear transition density as 



P\ O) 



P n x (r)±p p x (r). 



(5) 



On the other hand, pi (r) can be obtained using the 
same Bohr-Mottelson method, by deforming the IS and 
IV parts of the g.s. density 



P\ 



0(1) M - A , " LFg 
(r) — — (i) 
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dr 



(6) 



It is straightforward to derive in this compact approach [5] 
a consistent one-to-one correspondence between 5 p ( n ) an d 
So 



(i) 



which can be used to determine <Wi) values from 
values given by the DWBA analysis. If one assumes 



that the excitation is purely isoscalar and the neutron 
and proton densities have the same radial shape, scaled 
by the ratio N/Z, then S n = S p = So = Si. Therefore, 
any significant difference between So and Si would directly 
indicate a different isospin distribution in the structure of 
the nuclear excitation under study. 

In the present work, we have studied the elastic and 
inelastic 30 ' 32 S+p scattering data at 53 MeV/u [5] and 
34,36,38,40 S+p data [g ]|]fi]|n] ir2] at energies of 28 to 39 
MeV/u. The IS and IV contributions of the inelastic FF 
were considered explicitly to find out a systematic behav- 
ior of S\ along the Sulfur isotopic chain, passing by the 
magic number N = 20. Then, the folding + DWBA anal- 
ysis of the elastic and inelastic 18,20 O+p data at 43 MeV/u 
[13] and 22 0+p data at 46.6 MeV/u T4] has been done to 
find out the A-dependence of 6± in the Oxygen case. 

To have the accurate "distorted" waves for the DWBA 
calculation of inelastic scattering, the optical model (OM) 
analysis of the elastic data was done using the real folded 
potential [7] obtained with the density- and isospin de- 
pendent CDM3Y6 interaction [15] and nuclear g.s. densi- 
ties given by the Hartree-Fock-Bogoliubov (HFB) calcu- 
lation [16] . The imaginary part of the OP was taken in 
the Woods-Saxon (WS) form from the global systematics 
CH89 [17] . All the considered elastic data were well repro- 
duced (see Figs. [2] and with the depth of the WS imag- 
inary potential slightly adjusted by the OM fit (keeping 
the radius and diffuseness unchanged). The experimental 
reduced electric transition rate B(E2 f) was used in each 
case to fix S p value in the expression ^ for p\(r) which is 
further used in the folding calculation [5,7! ■ As the only 
parameter, S n was adjusted iteratively in the folding + 
DWBA calculation to fit the measured inelastic cross sec- 
tion. Since the CDM3Y6 interaction is real, the imaginary 
nuclear FF was obtained by deforming the imaginary part 
of the OP with So and Si values at each iteration step of 
the folding + DWBA fit to the inelastic data. In each 
case, the final set of deformation lengths So and Si was 
fixed only after the best-fit S n has been obtained. 

An earlier folding + DWBA analysis [7] of the same in- 
elastic 30-40 S-|-p scattering data, using inelastic FF given 
by the microscopic transition densities obtained in the 
Quasiparticle Random Phase Approximation (QRPA) [8], 
has shown that the neutron and proton contributions to 
the 2^ excitation in 3 °. 32 - 34 S follow approximately the 
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Fig. 1. The measured B(E2 f) transition strength (upper 
part), and the extracted isoscalar and isovector deformation 
lengths (lower part) of 2f states in Sulfur isotopes. The lines 
are to guide the eye. 



isoscalar rule [5] which implies So ~ 5\. The present fold- 
ing model analysis using the collective-model transition 
densities (0]) has shown about the same results (see So 
and Si values extracted for 3 °. 32 . 34 S in Fig. [I]). With the 
neutron shell becomes closed at N = 20, a significant 
"damping" of the neutron transition strength occurs and 
suppresses strongly the IV deformation length Si of the 
2f state of 36 S. In fact, Si is reaching its minimum as the 
neutron number TV approaches the magic number 20. The 
A-dependence of the IV deformation length Si is well cor- 
related with the A-dependence of the reduced transition 
rate B(E2 f) which also reaches its minimum at N = 20 
(see Fig. [IJ. The shell closure effect is so strong in this 
case that the proton and neutron QRPA transition den- 
sities needed to be scaled down by a factor of 0.63 and 
0.88, respectively, for a correct description of the mea- 
sured 2 + cross section in our earlier DWBA analysis [7] 
of the inelastic 36 S+p scattering. The contribution by the 
two valence neutrons in 38 S to the 2f excitation is quite 
strong and 5 n turned out to be larger than S p by around 
30% (see Fig. This difference between the proton and 
neutron transition strengths results on the IV deformation 
length Si larger than So by about 64% (see Fig. []}. Note 
that such a strong core polarization by the two valence 
neutrons in the 38 S case could not be fully accounted for 
by the QRPA calculation and the QRPA neutron tran- 
sition density has been scaled by a factor of 1.25 for a 
good agreement of calculated DWBA cross section with 
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Fig. 2. Elastic and inelastic 38 S+p scattering data at 39 
MeV/u [11] in comparison with the DWBA cross sections given 
by the folded FF. The dashed curve is inelastic cross section 
obtained in the isoscalar limit (S n = 6 P = So = 5i). 
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Fig. 3. The measured B(E2 f) transition strength (upper 
part) and the extracted IS and IV deformation lengths (lower 
part) of 2+ states in the Oxygen isotopes. The values for 24 
are deduced from the QRPA prediction by Khan et al. [19] . 
The lines are to guide the eye. 



the data [7J. Since the inelastic 40 S+p scattering data at 
30 MeV/u [12] consist of only two data points with signif- 
icant error bars, the extracted <Wi) values for 40 S might 
not be as reliable as those extracted for 30_38 S. The new 
measurement is, therefore, needed for an accurate estimate 
of the IV deformation length of 2^ state of 40 S. 

The present analysis of the inelastic 18 > 20 > 22 O+p scat- 
tering data has shown an interesting ./V-dependence of 5\. 
Except for the double-magic (N = Z) 16 0, where one 
has exact relation So = Si, the results obtained for 18 O 
already indicate a rather strong IV deformation for the 
2^ state of this nucleus (see Fig. [3]). We recall that the 
IS and IV deformation parameters of the 2^~ state of 18 O 
have been determined long ago by Grabmayr et al. [4] in 
a simultaneous DWBA analysis of the (p,p') and (n, n') 
inelastic scattering data at 24 MeV using the collective 
form factor ([2]). It is easy to deduce from the results of 
Ref. [4] the corresponding deformation lengths Sq w 1.1 
and Si m 2.6 ± 1.3 fm for the 2+ state of 18 0. These values 
agree reasonably with the results of our folding + DWBA 
analysis (see Fig. [3]) of the inelastic 18 0+p scattering data 
at 43 MeV/u [13]. Note that the 5i value extracted from 
our analysis has a much smaller uncertainty compared to 
that deduced from the results of Ref. [4] . The error bars for 
<5 (i) plotted in Figs. Q] and [3] were accumulated from the 
experimental uncertainties of the measured B(E2) values 
and ija,p') cross section. 



As already mentioned above, a double measurement of 
(p,p') and (n, n') inelastic scattering is not possible with 
unstable Oxygen isotopes, and our folding method is the 
only alternative way to determine 5%. The present anal- 
ysis (with a more consistent treatment of the IV part of 
the imaginary FF) has confirmed again the large IV de- 
formation length of the 2^ state of 20 O found earlier in 
Ref. 5j. With the IV deformation about three times the 
IS deformation (see Fig. [3]), the contribution by the Lane 
form factor Fi to the 2f excitation of 20 O amounts up to 
40-50% of the total inelastic cross section [5]. If we con- 
sider 20 O as consisting of the 16 O (or 18 O) core and four 
(or two) valence neutrons, then a large value of IV defor- 
mation length Si indicates a strong core polarization by 
the valence neutrons in the 2+ state of 20 O. 

In such a "core + valence neutrons" picture, it is natu- 
ral to expect that the l\ state of 22 O should be more col- 
lective and have a larger IV deformation length due to the 
contribution of two more valence neutrons. However, the 
inelastic 22 0+p scattering data at 46.6 MeV/u measured 
recently at GANIL [14] show clearly the opposite effect, 
with the (p,p') cross section about 3 to 4 times smaller 
than that measured for the 2+ state of 20 O at 43 MeV [13] 
over a wide angular range. The folding + DWBA analysis 
P3] of these data using the QRPA transition densities for 
the 2^ state of 22 O has pointed to a much weaker neutron 
transition strength compared to that of the 2+ state of 
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20 O. Given a significantly higher excitation energy of this 
state (1.5 MeV higher than that of the l\ state of 20 O), 
the newly measured inelastic 22 0+p scattering data were 
suggested [14) as an important evidence for the neutron 
shell closure at N — 14 or 16. 
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Fig. 4. The same as Fig. [2] but for the elastic and inelastic 
22 0+p scattering data at 46.6 MeV/u [14]. 



With the proton deformation length S v fixed by the 
adopted B(E2 f) value for the 2+ state of 22 [IS], the 
present analysis has given the upper limit for the neutron 
deformation length as 8 n < 1.2 S p (see Fig. [4]), which im- 
plies 5i < 1.1 fm. As a result, the extracted Si value is 
just slightly larger than So and the iV-dependence shown 
in Fig. [4] indicates that Si reaches its minimum at either 
N — 14 or 16. Based on the results obtained above for 
the Sulfur isotopes, we conclude that the deduced iV- 
dependence of Si for the 2^ states of Oxygen isotopes 
suggests the neutron shell closure at cither N = 14 or 
16. We note further that the B(E2 |) value predicted 
by the QRPA calculation [19 for the 2^ state of 24 is 
even smaller than that adopted for the 2j state of 22 O, 
a fact which could favor the shell closure at N = 16. We 
have also deduced Sq and eh of the 2^~ state of 24 from 
these QRPA results and they are nearly the same as those 
deduced from the (p,p r ) data in 22 case. If we take into 
consideration 4 MeV gap between the 2s!/ 2 and ld 3 / 2 sub- 
shells predicted recently by a consistent HFB calculation 
[2"0] . then it is also likely that the neutron shell closure 
occurs at N = 16. In any case, more measurements for 



24 O are highly desirable for a definitive conclusion on a 
new magic number N — 16 in the neutron rich nuclei. 

In summary, we have shown that the behavior of the 
dynamic isovector deformation of the 2^~ states in neutron 
rich nuclei is closely correlated with the evolution of the 
valence neutron shell. This interesting result emphasizes 
again the importance of (p,p') reactions measured with 
unstable nuclei in the inverse kinematics. 
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